Domain Structures in Fourth-Order Phase and 
Ginzburg-Landau Equations 



David Raitt 4 " 

Hermann Riecke 
Department of Engineering Sciences and Applied Mathematics 
Northwestern University 
Evanston, IL 60208, USA 



on ■ 

on : 
m ■ 



> 

m 
O 
O 
(N 
O 
^1" 
ON 

o 



o3 
Oh. 



13 



In pattern-forming systems, competition between patterns 
with different wave numbers can lead to domain structures, 
which consist of regions with differing wave numbers sepa- 
rated by domain walls. For domain structures well above 
threshold we employ the appropriate phase equation and ob- 
tain detailed qualitative agreement with recent experiments. 
Close to threshold a fourth-order Ginzburg-Landau equation 
is used which describes a steady bifurcation in systems with 
two competing critical wave numbers. The existence and sta- 
bility regime of domain structures is found to be very intri- 
cate due to interactions with other modes. In contrast to 
the phase equation the Ginzburg-Landau equation allows a 
spatially oscillatory interaction of the domain walls. Thus, 
close to threshold domain structures need not undergo the 
coarsening dynamics found in the phase equation far above 
threshold, and can be stable even without phase conserva- 
tion. We study their regime of stability as a function of their 
(quantized) length. Domain structures are related to zig-zags 
in two-dimensional systems. The latter are therefore expected 
to be stable only when quenched far enough beyond the zig- 
zag instability. 
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I. INTRODUCTION 

The formation of steady spatial structures in systems 
far from equilibrium has been studied in great detail over 
the past years, the classical examples being Rayleigh- 
Benard convection and Taylor vortex flow |l],^j . In quasi- 
one-dimensional geometries they usually share a com- 
mon feature: the stable structures that arise after the 
decay of transients are strictly periodic in space and if 
they are weakly perturbed they relax diffusively back 
to the periodic state. This relaxation has been inves- 
tigated experimentally in various systems ||,f|]|] and the 
results agree with theoretical results based on the phase- 
diffusion equation Jy|fl||l- This equation describes the 
slow dynamics of the local phase, the gradient of which 
is the local wave number. The structures are generally 



stable over a range of wave numbers which is limited (at 
least) by the Eckhaus instability The Eckhaus insta- 
bility is characterized by a vanishing of the phase diffu- 
sion coefficient and leads usually to a new stable periodic 
state with a different wave number. This process involves 
phase slips, i.e. the destruction (or creation) of one or 
more roll pairs in convection, say. The Eckhaus instabil- 
ity has been investigated in great detail [pil|p]JT2]Jl^| , in 
particular in Taylor vortex flow • 

Recently, it has been pointed out in theoretical work 
that a vanishing of the diffusion coefficient does not nec- 
essarily invoke a phase slip |p^ , [l7| , ^8| , pj|j20| . Instead, un- 
der certain conditions, the structure can become inhomo- 
geneous and evolve to a stable domain structure which 
exhibits different wave numbers in different parts of the 
system and which does not relax to a periodic structure. 
Such a structure is illustrated in fig{I]. The thick line de- 
notes the local wave number and the thin line a typical 
physical quantity such as the mid-plane temperature in 
convection. 

Experimentally, inhomogeneous structures have been 
found in a variety of systems. Most of them involve 
at least one time-dependent structure, e.g. counterprop- 
agating traveling waves (or spirals) pl| , ^2| , (turbulent) 
twist vortices amidst regular Taylor vortices p3|] , local- 
ized traveling-wave pulses in binary-mixture convection 
pij and steady Turing patterns amidst chemical traveling 
waves p5[ . Recently, however, domain structures involv- 
ing only steady convection rolls of two different sizes have 
been observed in Rayleigh-Benard convection in a very 
narrow channel ^(|. It has been suggested that these 
structures may be related to the phase-diffusion mecha- 
nism discussed above |i"6fl . 

Similar phenomena are relevant in two-dimensional 
systems. There, roll (or stripe) patterns usually undergo 
an instability at small wave numbers in which the phase- 
diffusion coefficient for perturbations along the rolls goes 
through zero. This leads to an undulatory deformation 
of the rolls which can grow into a zig-zag pattern. Such 
patterns are related to domain structures with the 'zigs' 
corresponding to domains of one wave number and the 
'zags' corresponding to domains of a different wave num- 
ber. 

Here we investigate steady domain structures from two 
points of view. For structures well above threshold a 
description of domain structures is obtained using the 
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jropriate 4th-order phase equation introduced earlier 
In sec.H w e perform a detailed study of this 
equation. We emphasize the bifurcation that leads to 
the domain structure and its connection to the Eckhaus 
instability. These results are compared with recent ex- 
periments p6| . 

Domain structures can also form close to threshold if 
the basic state becomes unstable to two periodic patterns 
with different wave numbers. If these two patterns arise 
at similar values of the control parameter, the result- 
ing competition can be described by a Ginzburg-Landau 
equation. In sec. Ill we introduce and investigate the ap- 
propriate equation. We present a detailed analysis of the 
stability regime of domain structures for periodic bound- 
ary conditions. In part of the parameter regime they can 
again be described by the above phase equation. Sur- 
prisingly, we find no trace of the breakdown of the phase 
equation in the stability behavior of the domains. In fact, 
the domains can be stable even if the local wave number 
is well below the neutral curve, where the phase equation 
clearly does not apply. 

Within the phase equation, domain structures are only 
stable if the total phase, i.e. the total number of cells, is 
conserved. In sec. IV we therefore study the Ginzburg- 



Landau equation with a spatially ramped control param- 
eter which allows the total phase to change. We find 
stable domain structures even in this general case and 
explore their regime of existence. In sec.|v| we briefly ad- 
dress the relevance of our results to two-dimensional pat- 
terns and discuss the zig-zag instability and the stability 
of the resulting zig-zag patterns in isotropic systems, e.g . 
Rayleigh-Benard convection and Turing structures p8| , 
and in anisotropic systems like electroconvection in ne- 
matic liquid crystals P9|, 



II. DOMAIN STRUCTURES WITHIN THE 
PHASE EQUATION 

In this section we consider steady one-dimensional pe- 
riodic structures away from threshold. The slow dynam- 
ics of long-wavelength perturbations of such structures 
can be described by an equation for the local phase of 
the structure. To derive the phase equation one assumes 
the structure to be close to periodic and expands the 
relevant fields, e.g. the fluid velocity, according to 

v(x,t) = v (e~V(^,r)) + ev 1 (e~ 1 ip{X,T)) + h.o.t., 

(1) 

where the slow variables X — ex and T = e 2 t, e <C 1 have 
been introduced Q. Inserting this expansion into the 
basic equations, e.g. the Navier-Stokes equation, yields 
at O(e) the solvability condition 



d T ip = V(q,R)d x tp. 



(2) 



Note that the phase equation (g) is nonlinear through 
the dependence of the diffusion coefficient on the local 



wave number q = dxf- In addition, we have indicated 
the dependence of P on a control parameter, such as 
the Rayleigh number R in convection. As long as the 
diffusion coefficient V(q, R) is positive, perturbations de- 
cay diffusively and the structure is stable with regard to 
long- wavelength perturbations. 

If T>(q, R) becomes negative the structure loses stabil- 
ity and perturbations in the local wave number grow. 
Usually this leads to a phase slip and eventually to a pe- 
riodic state with a different wave number. If, however, 
T>(q, R) becomes negative in the center of the band of sta- 
ble wave numbers, the perturbations can saturate and a 
domain structure arises. This is quite naturally the case 
when the neutral curve has two minima (see sees. Ill, IV 
below) [ p~7|j27| . It can, however, even occur if the neutral 
curve has a single minimum Jl7| , p^[ . 

To ensure well-posedness of the phase equation for 
T>(q, R) < it is necessary to retain higher order deriva- 
tives. Introducing an additional slow time r = eH and 
expanding the phase around a spatially periodic struc- 
ture with wave number q^i l p{^t t ) — QoX + e(f>{X,r), 
one obtains 



8 T cp =(D + Edx4> + F{d x ^) 2 )d 2 x 4> - Gd x( f>, (3) 



where 



e 2 D = V(q ,R), eE = d q V(q 0l R), 
F = d 2 V{q ,R), G>0. 



(4) 



In anticipation of the situation relevant to this paper, we 
assume D, E, and F to be 0(1). These coefficients can 
also be determined in a linear stability analysis of the 
periodic state from the behavior of the growth rate at 
small values of the Floquet parameter (cf. (^p| ) below). 
The wave number Q — dx4> = {q — <Zo)/ e satisfies 



d T Q = d 2 x (DQ + |q 2 + |q 3 - Gd 2 x Q). 



(5) 



This equation is the Ginzburg-Landau equation for a con- 



served order parameter Q |3 
are most easily understood 



The static solutions to (|5|) 
Dy rewriting the right-hand 



side as the equation for a particle in a potential, 



a'xQ = -o Q u{Q), 



Q^ + CQ 



(6) 



with the integration constant C. Usually T>(q) is con- 
cave downward implying F < 0. This yields a potential 
as sketched in fig.^a, which allows two kinds of solutions. 
The three extrema correspond to spatially periodic struc- 
tures, of which only the maximum yields a stable solution 
within the full equation (^). In addition, there exists a 
family of unstable, oscillating (in X) solutions, which 
correspond to structures with spatially modulated wave 
numbers. It contains, as a limiting homoclinic 
orbit which describes a system with a localized 'dip' (or 
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'hump') in the wave number and an associated ampli- 
tude depression @. By contrast, if F > (figj|b) two 
maxima are possible allowing also a heteroclinic orbit 
connecting two stable periodic structures. In this paper 
we concentrate on such solutions. 

The connection between the different solutions is re- 
vealed by a bifurcation analysis close to the onset of the 
instability of the periodic structure. Since can be 
solved exactly in terms of elliptic integrals it is possible 
to perform the expansion in the exact solution. This ap- 
proach emphasizes the importance of phase conservation 
and is demonstrated in |32| . 

It is easier, however, to expand the phase in (0) 



nAif* + rf(A 2 e irx + B 2 e 2 ^ x ) + h.o.t. + c.c 



(7) 



and derive an evolution equation for A\ . Here P 
n integer, and the amplitudes depend on the slow time 
T = r/H, 77 <C 1. No term independent of X arises in (Q) 
because (|3j) is invariant under (j> — ► cj> + const. At fifth 
order one obtains 

dfA, = -D 4 P*A 1 + P 2 {^± _ F.P^A^A, 



3G 



5 E, 



216 G 3 
where the coefficients 



l^i| 4 ^i, 



(8) 



D = D V + V 4 D±, E = E V + V 2 E 2 , F = F v +t 1 2 F 2 (9) 

have been expanded around the tricrital point 

D v = -GP 2 , El = 6GF V P 2 . (10) 

Thus, the instability, which is signified by a sign change of 
the diffusion coefficient, need not always induce a phase 
slip, as it usually would in the Eckhaus instability. In- 
stead, if T>(q) is negative in only a (small) region around a 
minimum, i.e. if F = ^d 2 T>(q) is sufficiently positive, the 
bifurcation becomes supercritical and the domain struc- 
tures, which correspond to an oscillatory mode of the 
'particle in the potential', are stable. In the subcritical 
case the domain structures persist up to the saddle-node 
given by 



D A 



6GP 



[E V E 2 - 3F 2 GP 



2\2 



(11) 



To make contact with experiments [ p6| , we model the 
diffusion coefficient T>(q) in the vicinity of a minimum in 



V{R, q) = V{R mi q m ) + d R V ■ {R 
+ \d 2 q V-{q-q m ) 2 . 



Rm) 



(12) 



At a Rayleigh number R, the parameters D, E and F, 
which are calculated by expanding T>(q, R) around the 
periodic solution at wave number q$, are then given by 



D = d R V ■ (R - Rm) + -d 2 V ■ (q - q, 
E = d 2 V ■ (q - q m ), 



d 2 V. 



(13) 



Here we have chosen R m such that T>(q m ,R m ) — 0. A 
typical phase diagram for such a case is sketched in figj|a. 
The thin line denotes the onset of the instability, whereas 
the thick lines give the locus of the saddle nodes. When 
raising R the instability first appears at q m and is su- 
percritical. Just before the transition the phase diffusion 
coefficient is still positive but it is very small over some 
range of wave numbers, thus irregularities in the struc- 
ture are expected to decay very slowly. For larger R and 
consequently away from the minimum the bifurcation be- 



comes subcritical. The width Sq = 2Py6G/d 2 V of the 

supercritical regime depends on the length L = 2tt/P of 
the system and vanishes as L — > oo. In the subcritical 
regime, there exist two different domain structures in the 
region between the Eckhaus instability and the saddle- 
node line, one of which is unstable. This is related to 
the result of Brand and Deissler who find two different 
solutions of this type numerically Q ■ 

In convection experiments in a very narrow channel 
Hegseth et al. have found stable domain structures with 
different wave numbers [^6| . Brand and Deissler |l6| sug- 
gested that they are related to the solutions discussed 
above, a position that is supported by the present re- 
sults. When raising the Rayleigh number Hegseth et al. 
find a transition to domains, which can be supercritical or 
subcritical, depending on the average wave number q m . 
Their results are shown in fig.^jb for comparison with the 
predictions by the phase equation. This figure plots the 
maximum and minimum wave numbers of the pattern as 
the Rayleigh number is increased for two patterns which 
differ in their initial, homogeneous wave numbers. If the 
initial wave number is close to q m there is a smooth (su- 
percritical) transition to a domain structure as indicated 
by the circles. On the other hand, if the initial wave num- 
ber is far from q m , the periodic pattern persists to much 
higher values of the Rayleigh number and then jumps to 
the domain structure. The data for this case is shown by 
the squares and is indicative of a subcritical bifurcation. 

Near q m and R m , our analysis suggests that phase dif- 
fusion should become very slow. This agrees with the 
experimental finding that the domain structures do not 
relax to a periodic pattern within the observational time 
scale when R is decreased just below R m J|^j3^] . From 
the above discussion the diffusion coefficient T>(q) is ex- 
pected to have a bimodal shape in q and a minimum at 
(i? m ,<7 m ), and to go through zero linearly in R there. 
This could be tested experimentally with methods simi- 
lar to those used by Wu and Andereck in Taylor vortex 
flow @. 

Within the framework of the phase equation the in- 
teraction between domain walls is purely attractive |3l| . 
Therefore, in regimes in which the phase equation applies 
domain walls annihilate each other as long as this is com- 
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patible with phase conservation. Thus, if the boundary 
conditions enforce phase conservation, the system will 
eventually reach a state with only a single domain-wall 
pair. If the boundary conditions do not conserve phase, 
e.g. in the presence of subcritical ramping in the con- 
trol parameter R |?3|,§, domain structures exist only as 
transients jl6|,[l9). In general, however, the interaction 
may also be repulsive, allowing domain-wall pairs to be 
stable without phase conservation. Such a situ atio n can 
arise near threshold. It is examined in section IV using 



the Ginzburg-Landau equation which is introduced in the 
next section. 



III. DOMAIN STRUCTURES WITHIN A 
FOURTH-ORDER GINZBURG-LANDAU 
EQUATION 

Close to threshold domain structures can arise if the 
basic state becomes unstable with respect to two pat- 
terns with different wave numbers, i.e. if the neutral curve 
for the appearance of a pattern has two minima. If the 
wave numbers are close to each other, such a system can 
be described by a single Ginzburg-Landau equation. It 
should be emphasized that in addition the minima must 
be due to the same mode. This is the case, for instance, 
in parametrically driven standing waves |l7],^7|]. Such a 
Ginzburg-Landau equation is fourth-order in space, re- 
flecting the two minima of the neutral curve, 



d f A 



D 2 d%A 



iD 3 d%A 



bid%A 



x 



tA-t\A\ 2 A. 



(14) 



The convective amplitude A gives a typical quantity, e.g. 
the vertical velocity in convection, via 

v s {x,z,l) = Se l ^A(X,f)f(2) + h.o.t.+c.c. (15) 

where X — 8x and T — S t are slow variables, <5 1, and 
f(z) is the appropriate vertical eigcnfunction. In order 
for ( |I4| ) to be asymptotically valid D 2 d 2 ^A and iD 3 d^A 

have to be of the same order as D^d%A. Periodic solu- 
tions of this Ginzburg-Landau equation and their stabil- 
ity have been studied by Proctor Q] . 

In the following we solve ( |T4| ) numerically for a finite 
system of length L with periodic boundary conditions. 
We are interested in particular in the behavior of the do- 
mains when the wave-number gradients become steep. To 
achieve this we scan across the full instability regime of 
the periodic pattern by changing the average wave num- 
ber of the domain structure. To this end we investigate 
the dependence of the system on its length L. To keep 
L, the nondimensionalized system length, as a parameter 
we scale lengths using the critical wave number q c rather 
than L, 

X = X/q c , q = q c (l + SQ), L = L/q c , 



A = Aql^Jb i /f 1 f = T/{D i qt 
D 2 = D 2 D 4 q 2 c , D 3 =D 3 D 4 q c , t = SD 4 qi (16) 

With this scaling the Ginzburg-Landau equation be- 
comes 

d T A = D 2 d 2 x A + iD 3 d x A - d x A + 5L4 - \A\ 2 A. (17) 

Either D 2 or £ could still be scaled away, leaving the 
quantity D 2 /Y, as the primary control parameter. We 
have chosen to retain both parameters to simplify scan- 
ning each of them across zero. As pointed out by Tucker- 
mann and Barkley, periodic boundary conditions for the 
physical quantities like v~ z imply non-periodic boundary 
conditions for the amplitude A J35[, 



A{X + L)e iL/s = A(X) 



(18) 



Thus the conserved slow phase is the total slow phase 
A$ of the pattern, not only that of the slowly varying 
amplitude and is given by 



A$ = L 



5 I Qdx. 
h) 



(19) 



For small 5, changing the length has a very strong effect 
on the amplitude A, since small changes in the physi- 
cal wave number q imply large changes in the reduced 
wave number Q. In particular, with these boundary con- 
ditions the reduced wave number can be scanned across 
by changing the length of the system. With periodic 
boundary conditions on A, for which Q = ^j- ■ n, this 
would not be possible. 

The neutral stability curve and Eckhaus boundary for 
( |l7| ) are shown in fig.[| for D 3 = in the limit of infi- 
nite system length. The Eckhaus boundary is obtained 
by a linear stability analysis of the periodic solution 
A = A\e %cix that yields the growth rate a of sideband 
disturbances with wave number Q ± P, 

<r = -P i -{D 2 + 6Q 2 )P 2 -Y; + Q 2 D 2 + Q i ± (20) 
v/(S - Q 2 D 2 - Q 4 ) 2 + 4P 2 Q 2 (D 2 + 2(P 2 + Q 2 ) 2 ). 

The Eckhaus curve is obtained from this in the limit 
P — > 0. In each of the two minima of the neutral curve 
the stability of the periodic pattern is governed by the 
usual Eckhaus criterion. Note that for Q = 0, i.e. at the 
local maximum of the neutral curve, the periodic pat- 
tern is unstable for all E. The width of this unstable 
region is AQ = y/-D 2 /6. The term "inner Eckhaus 
boundary" used below will refer to those parts of each 
Eckhaus boundary for which \Q\ < \Q m in\, while "outer 
Eckhaus boundary" will refer t o those par ts for which 
\Q\> \Qmin\, where Q min = ±y/-D 2 /2D i are the min- 
ima of the neutral-stability and Eckhaus curves. 

In a finite system there is no longer a continuum of al- 
lowed wave numbers, instead there is a discrete spectrum 
{Qn} of them. This also implies that P is quantized in a 
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similar manner. Using the chosen scaling (|l6|), P is con- 
strained to integer multiples of . In this case all allowed 
periodic solutions can be stable. Eq.(pO|) shows that for 
this to occur it is necessary that the wave numbers are 
spaced widely enough to straddle the unstable region in 
the center, i.e. D2/P 2 > —5/2, and sufficient that Q = 
itself is stable, which happens for D2/P 2 > — 1. The 
latter is the case if the adjacent allowed wave numbers 
Q = ±2ir/L lie outside the neutral curve. Again, a pe- 
riodic solution that is unstable at onset retains its insta- 
bility for all values of E. 

Keeping D2 fixed, we have as control parameters the 
total slow phase of the system A$, the system length 
L, and E, which can be considered a reduced Rayleigh 
number. The domains are characterized by the difference 
AQ between the largest and the smallest wave number 
in the solution which thus serves as an order parameter 
(cf. fig§. 

The previous work done on the phase equation (cf. 
sec.|n|) gives some indication of what to expect in these 
simulations. For large S and small D 2 the phase equa- 
tion , and therefore its amplitude equation (JsJ) , is valid. 
In that region of parameter space it is expected that a 
branch of unstable domain structures bifurcates from the 
inner Eckhaus boundary and undergoes a saddle-node 
bifurcation to a branch of stable domain structures (cf. 
figj^). On the other hand, for values of the parameters 
near the minimum of the neutral-stability curve, the solu- 
tions can be described using the phase equation appropri- 
ate to the typical Eckhaus case. Thus it is expected that 
there is an unstable branch that bifurcates subcritically 
from the inner Eckhaus boundary and combines with an 
unstable branch bifurcating subcritically from the outer 
Eckhaus boundary. In the following we concentrate on 
the transition region between these two types of behavior. 

To calculate domain structures and examine their sta- 
bility we solve the Ginzburg-Landau equation (^) with 
boundary conditions fll8| ) numerically. The solutions are 
found using Newton's method and their stability is deter- 
mined using inverse iteration. The average reduced wave 
number Q = L QdX/L is used as a control parameter. 
To change this it is most straightforward to alter L by 
changing the grid spacing dx (~ 0.25 in this section). The 
additional parameters are given by A$ w 113.7, 5 = 0.1, 
D2 = —1.0 and E chosen as indicated. We concentrate 
on the symmetric case D3 = 0. 

A set of four solutions at representative values of the 
two control parameters are shown in fig.js]. The aver- 
age reduced wave number Q remains constant in the 
columns while E remains constant in each row. The local 
wave number is marked by the heavy solid line, whereas 
Re(e lqaX A) with q c = 1 is denoted by the thin line. As 
is apparent in fig.[5]a,c the dominant effect of increasing 
the average wave number Q is a decrease in the length of 
the low-wave number region in the center. A striking fea- 
ture, which is particularly important in the next section 
is demonstrated in fig.^a,b; the local wave number is non- 



monotonic and exhibits spatial oscillations. These oscil- 
lations become increasingly pronounced with decreasing 
E. We will comment on the relevance of this below. 

A sequence of bifurcation diagrams displaying the tran- 
sition from the large-E to the small-E behavior is given 
in fig|| In figj^a E = 1.2. The branch of the stable 
domain structure is indicated by the bold line terminat- 
ing in the saddle-node bifurcation A at Q m 0.68. There 
it collides with an unstable branch of solutions that bi- 
furcates off the periodic solution (AQ = 0) at the inner 
Eckhaus boundary (cf. fig.| at Q w 0.43). The dotted 
line represents another branch of unstable solutions that 
arises from the outer Eckhaus boundary (Q « 1.12)and 
undergoes a phase slip as Q — > which changes A$ by 
2tt • S. 

Lowering the control parameter E to 1.0, introduces 
new features (fig.^b). The original saddle-node bifurca- 
tion A remains, but two new saddle-node bifurcations, B 
and C, have appeared in the stable branch. These new 
saddle nodes mark the ends of a loop that grows as E 
is reduced. In going through the loop, the wave num- 
ber changes from having a local maximum and two local 
minima in the low-wave-number domain (as in fig.^Jb) to 
having only a single minimum in that domain. Addition- 
ally, the branch from the outer Eckhaus boundary moves 
closer to the stable solution branch. 

Reducing the parameter to E = 0.9 (fig j^c) produces 
significant changes in the bifurcation diagram. The outer 
Eckhaus branch merges with the branch of stable domain 
solutions to produce two new saddle nodes, D and E. As 
E is lowered further, the saddle-node pairs C and D as 
well as A and E move together (fig.||d) and annihilate 
leaving only one surviving saddle node, B. 

Lowering E to 0.3 produces no qualitative change in 
the behavior of the system. There continues to be only 
a single saddle node, at the tip of the loop (fig.||e). An 
additional (dashed) line is indicated in this figure. It rep- 
resents another solution branch that undergoes a phase 
slip of 4-7T • 5 as Q — > 0. Comparing this analysis with 
that of the phase equation (^) it is tempting to associate 
this branch with a phase perturbation (f> of wave number 
P = ^-2 (cf. (0)). It is therefore expected to merge with 
the periodic solution for large Q. We were, however, un- 
able to continue it that far, since additional saddle-node 
bifurcations arise. In the wave-number region of interest, 
it merges with the branch of stable domain structures be- 
tween E = 0.3 and E = 0.25 producing two new saddle- 
node bifurcations, F and G, shown in fig.[j]f. The latter 
annihilates with saddle node B around E = 0.2. Upon 
lowering E further, additional saddle nodes form and an- 
nihilate pairwise. This is presumably due to the inter- 
action with additional branches corresponding to phase 
slips 2n ■ S ■ n, n > 2. 

The above numerical results yield the phase diagram 
shown in figj^. This picture shows one half of the pa- 
rameter space and is symmetrically continued for Q < 0. 
Stable domain structures exist in the region above and 
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to the left of the line formed by the solid circles. Each 
of them indicates a saddle-node bifurcation of the stable 
branch. Above S » 1.2 there exists only a single line of 
saddle-node bifurcations (A). The first line of new dots 
(B and C), to the left of line A denotes the loop first 
shown in fig||b. The two saddle nodes in the loop are so 
close together that they cannot be resolved in this figure. 
For smaller E (ss 1.0), a hump arises from the creation of 
saddle nodes D and E, at the merging of the stable branch 
with the 27r-phase-slip branch. These saddle nodes anni- 
hilate pairwise, leaving only one (B). At £ « 0.3 another 
hump appears, caused by the merging of the 47r-phasc- 
slip branch and the branch of stable domain structures, 
creating saddle nodes F and G. Again, two of the sad- 
dle nodes merge and annihilate leaving only one. Similar 
humps that appear in the saddle-node line suggest that 
this process continues as £ is reduced and Q approaches 
zero. 

The phase equation breaks down when the local wave 
number approaches the neutral curve. As fig.0 shows we 
have been able to find stable domain structures for val- 
ues of £ > —0.177. In this regime, the wave number is 
locally below the neutral curve. Thus the domain struc- 
tures exist and are stable well beyond the validity of the 
phase equation. Strikingly, this leaves no trace in the 
stability properties of the domain structure. It is, how- 
ever, responsible for the oscillatory behavior of the wave 
number (cf. fig.||). It is not clear whether or not the do- 
main structures exist all the way to the minimum of the 
neutral curve, occurring at £ = —0.25. 

For simplicity we have concentrated on the special case 
D3 = 1 36 . If D3 ^ then one of the two minima in the 
neutral curve is lower than the other. Thus, if the control 
parameter is raised adiabatically, domain structures will 
not arise spontaneously from the basic state as the critical 
parameter value is crossed. Instead a periodic pattern 
with a wave number corresponding to the lower minimum 
will arise. If, however, the control parameter is quenched 
into the region in which there are two separate stable 
regions, then domain structures are expected to form. As 
shown by Proctor 0] , the region of stable wave numbers 
that arises at the higher value of the control parameter 
is closed at the top rather than open as is the case here. 
We do not expect this to have a strong effect, however, in 
the regions of parameter space under investigation here. 



IV. DOMAIN STRUCTURES WITHOUT PHASE 
CONSERVATION 

Within the phase equation the dynamics of domain 
structures depends sensitively on the boundary condi- 
tions. In particular, the stability of a domain struc- 
ture consisting of a single domain-wall pair requires the 
boundary conditions to conserve the phase. It is therefore 
of interest to study which properties of such a domain 
structure are independent of the boundary conditions 



within the Ginzburg-Landau equation. This information 
will also give insight into the properties of domain struc- 
tures consisting of many domain walls. In such structures 
adjacent domain- wall pairs can annihilate even if the to- 
tal phase is conserved since nearby domains can grow or 
shrink to compensate for the loss of phase. 

To investigate the behavior of domains in such a situ- 
ation we use the fourth-order Ginzburg-Landau equation 
introduced in the previous section. Boundary conditions 
which do not conserve the phase can be obtained by in- 
troducing a smooth subcritical ramp in the control pa- 
rameter E at each end of the system 1 3^|| . Since the am- 
plitude at each end of the system goes then to zero, phase 
can be added or removed through the boundaries. Effec- 
tively, a ramp poses boundary conditions on the wave 
number |33|,||]. In this system it selects either of the two 
minima of the neutral curve lQ m j„. Since these wave 
numbers are consistent with the wave numbers found in- 
side each domain, the ramps mimic the effect of distant 
domain structures. 

To simulate the time evolution numerically we use a 
Crank-Nicholson scheme with dx sa 0.034 and nx — 3200. 
The weakness of the attraction between domain walls 
suggested such a small grid spacing to eliminate pinning 
by the lattice points. The initial conditions for all of the 
simulations performed in this section consist of a step 
function with three domains, two of high wave number 
surrounding one of low wave number. The initial condi- 
tions are varied by changing the width of the central (low 
wave number) domain. 

In order to confirm that this ramping is sufficient to 
eliminate phase conservation, we start by investigating 
large values of the control parameter, E = 100 (with 
D2 = —1.0). In this region of parameter space, the 
phase equation governs the dynamics of the pattern and 
the domain walls will have purely attractive interaction. 
As expected, all of the different initial conditions tested 
converged to a periodic pattern with wave number Q m in- 
This shows in addition that the numerical parameters se- 
lected are sufficient to eliminate pinning on the grid in 
this regime. 

We next consider the behavior of domains for lower val- 
ues of the control parameter E. These simulations were 
performed for E = 1.0 and D2 = —1.0. For these values 
of the control parameter, close to the neutral-stability 
curve, the dynamics is no longer expected to be governed 
by the phase equation. The behavior is indeed quite dif- 
ferent from that observed for high E. Here the system 
evolves to a discrete number of different states, as shown 
in fig.||, where the total phase in the system (between the 
ramps) is plotted versus time. As the initial width of the 
central (low wave number) domain is increased, the total 
phase of the initial condition is decreased, leading to the 
different lines shown in this plot. The upper lines are 
those with only an initially narrow region of low wave 
numbers. They converge to a final state with a single 
minimum in wave number. The next set of lines have a 
slightly wider region of low wave number and converge 
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to a final state with two minima in the wave number. 
Likewise the lowest set of lines converge to a state with 
three minima in the wave number. The corresponding fi- 
nal states are shown in fig|]. These results show that the 
width of the domain structures is quantized and the final 
states can be characterized by the number of minima in 
the wave number. 

We now turn to the transition region between the low- 
E regime, where domain structures exist, and the large- E 
regime, where they do not. In the next section, we will 
relate the stability regime of domain structures to that 
of zig-zag patterns in two-dimensional systems. For that 
application it is more natural to change \D2\ than E since 
the former is directly related to the roll-spacing of an ini- 
tial straight-roll pattern. These two control parameters 
are equivalent since ( p"7| ) has only a single control param- 
eter, the ratio D\jY>. We therefore study the dependence 
on Di instead of on E. 

In order to determine the value of the control parame- 
ter at which the various states disappear, we create initial 
conditions consisting of domain structures with different 
widths for D2 = —1.0 and E = 1.0. Then D2 is increased 
adiabatically and a new steady solution is determined 
using Newton's method. The values of D2 and the total 
phase A$ are recorded for each step. Figure [to] sum- 
marizes these results, plotting A$ against D2. As D2 is 
increased, the total phase in each state decreases. This is 
mainly because, the wave number in each of the domains, 
Qmin ; decreases in absolute value as | D2 | decreases while 
the length of the domains stays essentially fixed. Even- 
tually, each domain structure vanishes through a saddle- 
node bifurcation. These bifurcations occur at different 
values of the control parameter for different size domains. 
The value of D2 at which the bifurcation occurs increases 
for the domains with one to four minima and then de- 
creases again for the state with five minima. The investi- 
gation of larger domains was not feasible with our finite 
difference code. 

In the region of parameter space where the domain 
structures are stable, the interaction between the domain 
walls must change sign as a function of their distance. In 
principle, at large distances this interaction can be de- 
scribed using perturbation methods on heteroclinic or- 
bits |37 38 1 and in simpler systems it has been shown 
that a spatially oscillatory behavior of the domain wall 
can lead to a locking of adjacent walls |37]]. The asymp- 
totic behavior of the solution far away from the domain 
wall can be determined through a linearization around 
the corresponding periodic solution of (|l7|). For the pe- 
riodic solution with wave number Q m i n such an analysis 
reveals spatially oscillatory eigenmodes, but also modes 
that decay monotonically in space J32| . Correspondingly, 
in the numerical simulations the wave number is clearly a 
combination of oscillatory and monotonic decaying func- 
tions for small L^l- Thus, for adjacent domain walls to 
lock into each other the existence of an oscillatory mode 
is not sufficient; it has to dominate the monotonic mode. 

At D2 ~ —0.35 the spatial decay rates of the two 



modes are the same (for E = 1.0). For larger D2, the 
decay rate of the monotonic eigenmode is lower and it 
will dominate the oscillatory mode so the dynamics will 
be that predicted by the phase equation. For D2 < —0.35 
the oscillatory mode has a smaller decay rate than the 
monotonic mode so repulsive as well as attractive inter- 
action is possible. The oscillations are clearly illustrated 
in fig|5|a,c where it is apparent that raising E (lowering 
D2) makes the oscillations near the domain wall much 
stronger. This argument suggests that domain structures 
could exist up to D2 ~ —0.35 and that wide ones persist 
to smaller values of I-D2I than narrow ones. This would 
be in agreement with the results for the states with one 
to four minima, but it does not explain the behavior of 
the five-minima state. Thus, a complete understanding 
of the numerical result seems to require a detailed calcu- 
lation of the interaction. 

The locking described in this section does not occur 
within the fourth-order phase equation ( p|) despite the 
presence of the fourth derivative. As eq^q) shows the 
phase equation is effectively only second order when con- 
sidering steady solutions. If a sixth-order term Hd^Q 
were kept in the phase equation, oscillatory behavior of 
the wave number could be described by the phase equa- 
tion for sufficiently large H . For this situation to be 
asymptotically valid both D and G have to be vanish- 
ingly small at the same time, which is, however, not to 
be expected in a typical physical system. 



V. DOMAIN STRUCTURES IN 
TWO-DIMENSIONAL SYSTEMS: ZIG-ZAGS 

The results discussed in this paper can be extended 
to describe certain aspects of two-dimensional patterns. 
A common instability of two-dimensional roll patterns is 
the zig-zag instability. It leads to an undulatory deforma- 
tion of the pattern along the rolls. During the nonlinear 
evolution the undulations become steeper and domains 
of parallel rolls arise. They are separated by domain 
walls within which the orientation of the rolls changes 
rapidly. Thus, they can be characterized by a wave vec- 
tor (q x ,Qy) with q y positive and relatively constant in a 
'zig'-region and negative in a 'zag'-region. The behavior 
of q y is therefore very similar to that of Q in the domain 
structures. Single domain walls between zigs and zags 
have been studied previously p7| . Here we concentrate 
on pairs of domain walls, which form a domain structure, 
and the interaction between them. 

In isotropic systems, such as Rayleigh-Benard convec- 
tion, the zig-zag instability arises generically when the 
wave number becomes too small, since the zig-zag pat- 
tern leads to an increase in the local wave number. While 
these patterns have not been observed to saturate in 
Rayleigh-Benard convection |4q[ , numerical simulations 
| 28| suggest that they can be stable in chemical systems 

[39pl. 
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Zig-zag patterns are also of great importance in 
anisotropic systems such as electro-hydrodynamic con- 
vection in nematic liquid crystals. There the long 
molecules exhibit a preferred direction and thus define an 
axis of anisotropy of the fluid. In these systems there are 
two possible primary bifurcations to roll- like structures. 
The first is to a parallel-roll pattern, in which the rolls are 
oriented normal to the axis of anisotropy and have wave 
vector (5^,0). Such a pattern is therefore commonly re- 
ferred to as a "normal-roll" pattern. The second type of 
primary bifurcation is to rolls oriented at an angle to the 
director. Due to reflection symmetry such rolls can be 
oriented in two directions with wave number (q x ,±q v ). 
This allows the formation of zig-zag patterns p{| ■ 

In a simplifying approach, it is possible to describe zig- 
zag-patterns using thephase equation with Q corre- 
sponding to q y |tl 42 ,43]. Thus, only variations in q y are 
considered, while q x is taken to be constant. As discussed 
above, such a description contains in general no stable 
zig-zag structures; instead there is a coarsening of the 
pattern until there are no domain walls left and the pat- 
tern consists of only one type of rotated rolls. Only if the 
boundary conditions enforce phase conservation, there 
will be a single domain wall left (a pair of domain walls 
for periodic boundary conditions) which corresponds to 
a 'zig-zag'-structure ('zig-zag-zig'-structure for periodic 
boundary conditions). 

In two-dimensional experimental systems the total 
phase is usually not conserved. It is therefore again of 
interest to investigate the stability of zig-zags in the ab- 
sence of phase conservation. Thus, we consider the ap- 
propriate Ginzburg-Landau equation. For isotropic sys- 
tems (as in Rayleigh-Benard convection) it is given by 



d T A 



■(id. 



x 



d 2 ) 2 A 



XA- \A\ 2 A, 



(21) 



while for anisotropic systems (as in convection in nemat- 
ics) it reads 



d T A = (& 



■ x - iZd x d Y 



2 1 Wdy — dy)A + XA — \A\ 2 A. 



(22) 



In the anisotropic case, the Lifshitz point, at which the 
pattern at onset changes from normal to oblique rolls, is 
given by W = -ZQ X . 

Concentrating on solutions which are strictly periodic 
in X, A = Ax(Y, T)e lQ * x , one obtains eq.(p7|) with 



D 2 = 2Q X , £ = A - Ql for © (23) 
D 2 = W + ZQ X , £ = X-Ql for (f|). (24) 

Note that Y and Q y in the two-dimensional systems cor- 
respond to X and Q in eq. (|l7|) and that now periodic 
boundary conditions for the physical quantities imply 
periodic boundary condtions for A\. Thus, the domain 
structures discussed in sec. Ill are indeed one-dimensional 
analogs of zig-zag patterns. 



Zig-zag patterns have been studied in some detail for 
anisotropic systems 

mill],||,||. TheSe studies found 
that there is no continuous transition from normal rolls 

to zig-zags of small amplitude. For sufficiently large am- 
plitude, however, stable zig-zags arc found to exist in is- 
lands within parameter space. Within these islands, the 
zigs and zags always orient themselves in the direction of 
maximal growth rate. Presumably these islands are re- 
lated to the discrete set of domains with different widths 
discussed above. 

Our results suggest a mechanism to explain the exis- 
tence of these islands of stability. As eqs.( pl[ - |2^ ) show, 
these results apply to both the anisotropic and the 
isotropic case. The implications of the results shown in 
figjiol are presented in figs.|ll|a-c. There the loci of the 
saddle-node bifurcations of the different domain states 
are given in the A — Q^-plane (medium- weight lines). 
These are the lines with = constant, where the ap- 

propriate constant for each domain structure is obtained 
from fig.[lO| In addition the Eckhaus curve and the onset 
of the zig-zag instability are pictured. The figures differ 
in the values of Z and W , with figjll](b) correspond- 
ing to the isotropic case. In each figure, the behavior of 
the patterns is similar. Between the Eckhaus instability 
(thin dotted line) and the zig-zag instability (thin dashed 
line) straight-roll patterns are stable. As the zig-zag in- 
stability is crossed they become unstable to undulations. 
Close to the onset of the instability the coarsening dy- 
namics predicted by the phase equation will be observed. 
This results in a pattern that will eventually consist of 
oblique rolls in only one direction. Their wave number 
will be close to that with maximal growth rate and they 
will be stable. 

If, however, the pattern is quenched deeper into the 
zig-zag-unstable regime the resulting zig-zag pattern will 
be able to persist, as the oscillations in Q x can lock into 
one another. The analysis of the spatial decay rates sug- 
gests that it is necessary to quench at least beyond the 
thickest solid line, which corresponds to D 2 < —0.35. 
The numerical results indicate one has to quench all the 
way beyond the medium- weight solid line, which corre- 
sponds to the four-minima state to find stable zig-zags. 
The different medium-weight lines denote the stability 
limits of 'zig-zag-zig' structures with different lengths of 
the 'zag'-domain. Note that these lines merge with the 
line denoting the zig-zag-instability on the neutral curve 
(thin solid line). Except for W = the merging occurs 
therefore at wave numbers for which the straight rolls are 
unstable with respect to the Eckhaus instability. Con- 
sequently, for small A stable, locked domain structures 
arise only at wave numbers Q x for which the straight 
rolls are unstable. Although the Eckhaus-instability of 
straight rolls does not necessarily imply that of the zig- 
zag-structures, it points to the importance of instabilities 
which involve the X- as well as the Y-dependence and 
suggests that these may be of particular importance for 
small A. In our analysis there is no allowance for such 
two-dimensional instabilities. For some of the domain 
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structures they have been investigated in the anisotropic 
case |45|j4l| , ^ , ^3| . In these studies stable zig-zag struc- 
tures have been found in certain parameter regimes. 

A typical 'zig-zag-zig' solution is illustrated in figjl^a 
for D2 = — 1 and £ = 0. It takes the one-dimensional 
pattern A(X) shown in figjl^b and extends it into two 
dimensions by plotting Re(A(X)e l ® yY ) with constant 
Q y = 0.7. The lines represent the zeroes of the pat- 
tern, while the dark spots are maxima and minima. This 
pattern is stable. A signature of the oscillations in Q x is 
the appearance of undulatory deformations in the zig-zag 
structure near the domain walls. 
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VI. CONCLUSION 

We have considered two approaches to investigat- 
ing the existence and stability of domain structures 
in pattern- forming systems: a phase equation and a 
Ginzburg-Landau equation. Within the phase equation 
the structures typically arise through a subcritical bi- 
furcation. This approach yields qualitative agreement 
with various experimental results in slot convection [p6| . 
As a conclusive test we suggest a measurement of the 
wave-number dependence of the phase diffusion coeffi- 
cient T)(q), which is expected to have a bimodal shape. 
The phase-equation approach becomes invalid when the 
transition region between adjacent domains becomes too 
narrow, i.e. when the gradient in wave number becomes 
too steep. This is in particular the case near onset. 

In the second approach we used a Ginzburg-Landau 
equation to describe the competition between patterns 
with different wave numbers close to onset. Such a com- 
petition naturally arises in a system with a double- welled 
neutral-stability curve, as is, for instance, possible in 
parametrically driven standing waves ]l7|Jl8|^7t . This 
approach allowed the investigation of domain structures 
and their stability even for large wave-number gradients. 
We found domain structures are stable over their entire 
region of existence. Strikingly this region is very intricate 
due to the interaction with unstable structures possibly 
arising from the outer Eckhaus boundary. 

The domain structures exist well beyond the regime of 
validity of the phase equation. A characteristic feature 
of those which are not described by the phase equation 
is the oscillatory behavior of their wave number. This 
allows adjacent domain walls to lock into each other, 
forming domain structures which are stable even without 
phase conservation. We have numerically studied the de- 
pendence of their regime of stability as a function of the 
length of the domains. 

Zig-zag patterns can be considered as two-dimensional 
analogs of domain structures. Therefore, in systems in 
which the phase is not conserved, they are not stable at 
onset of the zig-zag instability. Only further in the zig- 
zag-unstable regime can the locking prevent the pairwise 
annihilation of zigs and zags. 
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FIG. 1. Typical domain structure solution. The thick line 
gives the local wave number, while the thin line represents a 
typical physical quantity like the mid-plane temperature. 

FIG. 2. Potentials U(q) for (a) F < (V(q) concave down) 
(b) F > (D(q) concave up). 

FIG. 3. Comparison of (a) phase diagram obtained from 
the phase equatio n (J8J ) with (b) experimental results obtained 
by Hegseth et al Pqj T 

FIG. 4. Neutral stability and Eckhaus curves for D2 = — 1 
and D4 = 1. 

FIG. 5. Four solutions obtained at different values of the 
control parameters Q and E. The local wave number is 
marked by the heavy solid line, whereas Re(e tqcX A) with 
q c — 1 is denoted by the thin line. The order parameter 
AQ is also indicated. 

FIG. 6. Bifurcation diagrams for various values of E: a) 
E = 1.2, b) E = 1.0, c) E = 0.9, d) E = 0.8, e) E = 0.3, 
f) E = 0.25. The branches of stable domain structures are 
denoted by thick lines. For details see text. 

FIG. 7. Phase diagram: domain structures exist and are 
stable above and to the left of the saddle-node line denoted 
by solid circles. The letters denote the loci of the respective 
saddle-node bifurcations shown in fig|| The dashed and thin 
lines give the neutral and the Eckhaus curves, respectively. 

FIG. 8. Temporal evolution of the total phase for different 
initial conditions in the absence of phase conservation. The 
corresponding final states are shown in fig^. 

FIG. 9. Coexisting stable final states obtained from the 
time evolution shown in fig.^. The thick line gives the local 
wave number in the bulk. Note that the amplitude goes to 
zero toward the boundaries due to the subcritical ramping. 



FIG. 10. Regime of existence of bound pairs of domains 
walls for E = 1. The calculation was performed with up 
to 25,600 points and an appropriate grid spacing down to 
dx w 0.016. The ramped part is up to 3,200 points wide on 
each side. 
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FIG. 11. Extension of one-dimensional results of fig. 10 to 
zig-zag patterns for different values of W and Z in eq.(|22j)7(a) 
W = -0.25, Z = 1.5, (b) W = 0.0, Z = 2.0 (isotropic case), 
and (c) W = 0.25, Z = 1.0. To the left of the medium-weight 
lines zig-zags of different widths are locked: four-minima state 
(solid), one-minimum state (dotted). At the thickest line 
the monotonic and the oscillatory mode have the same decay 
rate. Between it and the zig-zag instability (thin dashed line), 
there are no stable straight-roll or zig-zag solutions within 
this framework. Thin lines: neutral-stability curve (solid), 
Eckhaus instability (dotted). Between the zig-zag instability 
and the Eckhaus instability straight-roll patterns are stable. 

FIG. 12. Two-dimensional analog of domain structure: 
zig-zag pattern (a) and the domain structure (b) from which 
it was generated. Both the wave number (thick line) and 
Re{Ae lq " x ) with q c = 1 (thin line) are shown in (b). The 
parameter values are E = 0, D2 = —1, Q y = 0.7. 
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